


Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1980-06 


An interactive code for a pressurized water 
reactor incorporating temperature and Xenon feedback 


Heath, Gregory Garver 


Monterey, California. Naval Postgraduate School 
http://ndl.handle.net/10945/17616 


This publication is a work of the U.S. Government as defined in Title 17, United 
States Code, Section 101. Copyright protection is not available for this work in the 
United States. 


Downloaded from NPS Archive: Calhoun 


Calhoun is the Naval Postgraduate School's public access digital repository for 
(8 DUDLEY research materials and institutional publications created by the NPS community. 
«ist Spe Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS'‘s first 


INN KNOX appointed — and published -- scholarly author. 

| LIBRARY Dudley Knox Library / Naval Postgraduate School 

411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 


r x al | View 
b rs c\ 
+ ‘ | ' rT ny! 1 iY ie 3 
i re 
Ly vl i 
rT 4 _ 
1 ' Ls 
, 1 e 
i] LP es 
a . 
+ , Le 
Fe » 
J a Ll 






























































- % } a i] , 
pO i F a7 a 
f 
. ’ , O 4 a Cy re 
. fs ie aS a 
. ' os Ww - 
i 
rs 
e i ‘+ a ‘ a i 
“ 4 4 y r ree 
F * oe ta 
iT * 
. 4 ‘ i= fl Fae, 
s , le 
i ‘ Ye * 1 ‘i 3 
% cI Py as C5) ¥ ivi 
Fi > e 7 ‘an ¢ t ¥ af Pe: Fie ‘ Fi 
a » ‘ _ a , bed 
; _ SEAR LCST rN 3 2) \ ne 
H C H e °,. = r eo oa | 2 
. ’ , . Md i i ' ie 
: i ae Nar at ‘ oes ee 
. rd ‘| 
ra » a . Con es rt Sy MY ; 
1 ey . ee A “a it > bd ‘s 
. 5 E Py at ¥ te € ' sah oye iz oe — 
1 O * Par ee ! i ora ‘ 
. ] i el 
rs 2 ry bd be a - ‘ a t 
. a td ‘ u Ly a v m1 HT 
ze Ny “1 y ir 5 s 
sh * = din 
. . a i ae | re ie te 
® my A rN) 
a an 4 a . 
) Ls - 
4 4 A O b OT tm r Be 
J J ™ 
F bd > " (3 a = 
Sy Ly o [ nt 4 
Ci eos eo s 
; : = ae ha ' 
5 » 
| c cA aha. ct 
. . 4 R ry Y i) 
. Mae L iu . . ae J tod 
a = Ps a ¢ 
, -_ 
‘ bs ' 
. a 
, 4 
€ - ' 
x 
A a 
. 
4 s J e 
oe 
* 
i * a 
ro hI r 
ei 
r A Py 
ry 
7 
A 
RB = i 
Py - 
f 
ry rl 
. 
. D 
H 
ea Ld ° os 
« 
a 
ea 
. 
5 
‘ 
' 
‘ @ rd 
ri ba i} 
ro 
Fo 
* 
Lf . Pd 7 He 
a 
mo 
a 
5 ad 
V7 
’ F 
uJ @ i- 
PJ 
« 4 
P t) 
rd 
p 
ry r 
? 5 
Ld Ci al s 
‘ 
3 
ff 
f 4 Pa 
UO a 
<¢ : 4 fy ’ “; 
ba ad 
ty a 
Ps rd ra of 
= 
a PM a al Pak 
J a 
Pd Ps - - 
ny ' ry ty ra od a ad 
r ae * 
| | F / va) ld 
r Ps is ae al 
s a rT s — = 
12 ta 7 a a - 
5 te - bad s a hI i 
x ‘ vr FI 
a Uy Bi ™ aa 4 fof 
O , 5 el a) Py P LT ig] 
® q - ro f i eae ¢ 
e - ry = te 
ys i ble . oe a ue? ary 
. ? 5 it ° f > bd elas F ‘ 
; “ re Ce) A Ps = y i gf ten) 
Pv) r ao PY a —-— 
6 a Ce i ts < ar reek 
F % . 7 ‘ i ” -; i F a vise ety 
r we be ra “. « ° sa nm tire ae 
e a a ? iW u ry es ue = 
« as F re fie | x de *. os 
ry = 5 r 
r « me f ‘ nat « “™ ary 9 wpe 
f Fy «a , ° Pro ad ” hee oy 
A a} FS ia 2 ad = wf ei 
” a oe Per aT ot SL ad Me aera her Bn, 
o ide » - 
3 [ ce : = = . > 
; , Di : # % pL Fj aA i y eT ag Pf 
mo Pi ee ¥, a ae a 
. a | rie a oe FI re poor ia ae wr Ci 
aos o fv) ‘ y eg AE ig Lrered 
a) Peds er i ve deh: err * = Wy a of Bs 
; ts 4 4 we aeae “oe ee was Cy Seatac dain 
r wu CD eT id een? pried ang i ot, b 4" Sf 
4 § he ee an oF e oY] ¥ Pr fd a + ree neg Me 
« » Pe i oe | F Pn od, Le ey of J ie 
. a e 4 Prd a re ee a dod 
Pr Pe hh Gy bit ‘7 5 
Ly : en ra : 1 A JF ge 
bs a y aD ® » ,& oT ee] a ee i” MM: ated 
* ‘ yr f an 24, i ey ee ¥, 
' Cf ob fy yr Te) Lt TF Sb) eke ee ae wer Ld 
f : H et ir BR er y / D 1 . oe 
as ” oe 7 y =iyh a ait lle Melia hat 
‘ 7: i 
ri ° usae Ps ae » hy ee a Poe , i ms Me 
Fi n ; oh i "yi nee i 
Pa) ae 7 * « r Pt 4 Li “ « Un at TF ats ie ie ie ae ras 0% ba 
os i i . t ‘3 
we p Ps 7 vey arty yr? oa _ us hes ty $5} 
¢ ® s a) b U s 4 y A i if 
, ina C rT ¢., : | ’ ret a aT ad = qe a ey ee abe I 
. * ° , Ey in ao Po . ti 
U ' oe wr e n " ee ae ay Prez 5 em pul’ are aps Tee as 
4 Ne ag CS re ay 0 F “ 
p $ “os A S een ng od ana Ta | : ‘“ sy ed Un 
J . i 
a ay ry ed rw} ‘ 4 ¥ aa) vi ty ia ie a Ppt 
Pa n ' oe ae a r A . “a r i hs oe oe 
f : " Ls ry iy i ve i i A ge 
1, i lh aa a yar) Ny 1 re Hy ‘y yee A 
r ’ F ri . % U0 ey ls 
A - A, > me he \ or ie ASN 
; tow Ae Te Os “ a Pe ae DF 
¥ t P 1G fen eG 
F 1 PA ‘ Pe 4 115 eft +4 Pee: He. 
I J Pore.) a 
4 in ree f i ee hae Ty m Pee 
) { bY a Se * i i q ry Aye ay ae i BY A 0 


DUDLEY KNOX ‘BA, ARY 
NAVAL POSTGRADUATE SCHOOL 


MONTEREY. CALIF 23940 








fi 











NAVAL POSTGRADUATE SCHOOL 


Monterey, California 





vray 


AN INTERACTIVE CODE FOR A PRESSURIZED WATER REACTOR 
INCORPORATING TEMPERATURE AND XENON FEEDBACK 


by 


Gregory Garver Heath 


June 1980 


Thesis Advisor: P., J. Marto 





Approved for public release; distribution unlimited 


1196581 





Unclassified 


SECURITY CLASSIFICATION OF THIS PAGE (Wren Dota Entered) 


REPORT DOCUMENTATION PAGE 


2. GOVT ACCESSION NO 








READ INSTRUCTIONS 
BEFORE COMPLETING FORM 


3. RECIPIENT'S CATALOG NUMBER 


5. TYPE OF REPORT & PERIOD COVERED 


Master's Thesis: June 1980 






TITLE (and Subtitie) 


An Interactive Code for a Pressurized Water 
Reactor Incorporating Temperature and Xenon 
Feedback 


AU THOR/(a) 





4 




















6. PERFORMING ORG. REPORT NUMBER 








7. 6. CONTRACT OR GRANT NUMBER/(e) 






Gregory Garver Heath 






10. PROGRAM ELEMENT, PROJECT, Tasx 


9. PERFORMING ORGANIZATION NAME AND ADORESS 
AREA & WORK UMIT NUMBERS 


Naval Postgraduate School 
Monterey, California 93940 















12. REPORT OATE 


ine_ 1980 


13. MUMBER OF PAGES 
a3 


1 CONTROLLING OFFICE NAME ANDO ADDRESS 


Naval Postgraduate School 
Monterey, California 93940 


— MOMITORING AGENCY NamE &@ AOORESSE/( II diferent jram Controtling Office) 18. SECURITY CLASS. (0/ (hte report) 








Naval Postgraduate School Unclassified 
Monterey, California 93940 


$e. DECLASBSBSIFICATION/ DOWNGRADING 
SCMEOULE 


16. OISTRIBUTION STATEMENT (0/ thie Repert) 


Approved for public release; distribution unlimited 


17. BCISTRIGUTION STATEMENT (af (he anarcact entered in Bleek 20, if different trem Report) 


18. SUPPLEMENTARY NOTES 


19. KEY wWOROS (Continue on reveree side if nececoary and identify by block number) 


Pressurized Water Reactor, Xenon Feedback, Moderator Feedback 





20. ABSTRACT (Continue an reverse side if necescary and identify by bfeck mamber) 
An interactive computer model of a highly enriched pressurized water reactor 
was developed, using the applicable plant parameters from Shippingport Atomic 
Power Station. The point reactor kinetics equations for one delayed neutron 
precursor group were linearized using small perturbation theory. The model 
included both moderator and Xenon-135 reactivity feedback effects, as well as 
an automatic reactor protection and average reactor coolant temperature control 
system. The thermal response of the model plant was simulated for normal 
operating transients induced either by control rod or turbine load changes. The 










DD , ee 5 1473 EDITION OF 1 NOV 68 18 OBSOLETE Unclassified 





Uncilassiried 
Gecumly CL AsSiPIC ATION GF TIS PAGErwren Nore Entered: 








post-shutdown Xenon transient response was also modeled. The interactive 
program was coded in FORTRAN-IV language, and the simulation program was 
coded in IBM CSMP-III language. 


DD : 
| a ; Unclassified 


— B=) 2s Bae .hUl ee ES es ictal 





Approved for public release; distribution unlimited 


An Interactive Code for a 
Pressurized Water Reactor 
Incorporating 
Temperature and Xenon Feedback 


by 


Gregory Garver Heath 
Lieutenant Commander, United States Navy 
B.S.M.E., United States Naval Academy, 1971 


Submitted in partial fulfillment of the 
requirements for the degree of 


MAStER@ OF SCIENCE INSENGINEERING SCIENCE 


from the 


NAVAL POSTGRADUATE SCHOOL 
June, 1980 





OUOLEY KNOX 


iBFiAiR 
NAVAL POSTGRan: 

» ADU oie 
MONTEREY, Cais ots! 


ABSTRACT 


An interactive computer model of a highly enriched pressurized water 
reactor was developed, using the applicable plant parameters from the 
Shippingport Atomic Power Station. The point reactor kinetics equations 
for one delayed neutron precursor group were linearized using small 
perturbation theory. The model included both moderator and Xenon-135 
reactivity feedback effects, as well as an automatic reactor protection 
and average reactor coolant temperature control system. The thermal 
response of the model plant was simulated for normal operating transients 
induced either by control rod or turbine load changes. The post shutdown 
Xenon transient response was also modeled. The interactive program was 
coded in FORTRAN-IV language, and the simulation program was coded in 


IBM CSMP-III language. 
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I. INTRODUCTION 


Computers have been used in the design and analysis of nuclear 
reactors since the inception of reactor technology, and many codes have 
been developed. More recently, digital computer programs have been used 
as learning devices for nuclear engineering students. The formulation 
of the computational problem of predicting the behavior of a nuclear re- 
actor has been greatly facilitated with the use of high level computer 
languages. 

The purpose of this work was to develop a computer assisted learning 
device to be used by students taking nuclear engineering courses at the 
Naval Postgraduate School. The program graphically displays the simulated 
kinetic and thermal transient responses of a pressurized water reactor 
power plant. The reactivity feedback effects of Xenon-135 poisoning and 
moderator temperature are incorporated into the model. Normal operating 
transients, starting from a steady state critical condition, can be 
induced by ramp changes in either control rod position or turbine load. 
The initial reactor power level and reactivity change mechanism are 
chosen by the program user. For user convenience, these inputs are 
prompted and entered interactively, after which the simulation is run 
with no further user action required. 

A pressurized water reactor (PWR) was modeled as it is the most 
common type of reactor used for power plant application. As a realistic 
reference, the model design incorporated the applicable characteristics 
of the Shippingport Atomic Power Station, whose nuclear parameters were 
the most compatible, of available PWR core data, with the assumption of 
a highly enriched core made in the model design. 


9 





The model was developed by treating the reactor kinetics and other 
plant component thermodynamic relationships as transfer functions [1]. 
The point reactor kinetics equations were linearized using a Taylor ex- 
pansion or small perturbation technique [2]. This same perturbation 
technique and lumped parameter analysis were applied to the plant's 
linear differential heat transfer and Xenon-135 equations [3]. The 
individual transfer functions were then integrated into an overall plant 
block diagram. The program also features a reactor protection and 
average reactor coolant temperature control scheme. 

IBM CSMP-III language was chosen to formulate the model simulation 
as it has inherent routines to invert the transfer functions to the time 


domain. 
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II. MODEL DESIGN CONSIDERATIONS 


Lumped parameter analysis and first order perturbation theory were 
used throughout the model development. All the nuclear variables which 
were used in the reactor kinetics and Xenon-135 decay equations were 
considered to be averaged values of the variables over the neutron 
energy spectrum. Similarly, the thermodynamic variables which were 
used in the plant's heat transfer equations were considered to be 
averaged values over the volume of the particular component. 

Because of these approximations and other assumptions made in the 
model development, the simulation should not be considered a design or 
Stability analysis. Instead, the simulation shows the model's large 
scale reactor kinetic and thermal trends during normal operating 
transients. 

For similar reasons, although certain plant parameters of the 
Shippingport Atomic Power Station were incorporated into the model, the 
simulation cannot be considered to reflect the operating characteristics 
of this power plant. The Shippingport core is a seed and blanket type 
whereas the model core is a uniform mixture of highly enriched fuel 


and support materials. 


A. PRESSURIZED WATER REACTOR 

The pressurized water reactor (PWR) is the most widely used reactor 
type in central power plant applications and the only type currently 
used in naval propulsion. 

A simplified schematic of the modeled PWR plant is shown in Figure 


1. The modeled reactor contains a highly enriched core which is light 
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Figure 1. Schematic of the Model Plant 


water moderated and cooled. The modeled plant contains two closed loop 
thermodynamic systems coupled by a heat exchanger. 

In the primary system, heat generated from thermal fission is trans- 
ferred to the coolant as it passes through the reactor, raising the 
coolant's temperature. The high temperature coolant leaving the reactor 
is circulated by a pump through tubes inside of a heat exchanger where 
it gives up some heat to a secondary fluid. A pressurizer maintains 
the primary coolant at a sufficiently high pressure so that the bulk 
temperature of the coolant is kept below the saturation temperature. 

On the secondary side of the heat exchanger, the temperature of the 
entering feedwater is raised to the boiling point and saturated steam 
1S produced. The steam is delivered to a turbine, is condensed and 


returned to the heat exchanger, completing the second closed loop. 
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B. MODEL AND PROGRAM FLEXIBILITY 


The model was designed to simulate only normal operating transients 
from an initial steady state condition. Reactor accidents and startup 
were not considered in the model development. However, the model's re- 
actor control module will simulate either a full or constant insertion 
of control rods if certain parameters are exceeded. This feature was 
not incorporated for accident analysis, but to keep parameters within 
the limits of the assumptions made in the model development. 

The applicable Shippingport plant characteristics are fixed in the 
simulation program. User inputs are limited to choosing initial power 
level, from which other initial parameters are adjusted, and the plant 
perturbation mechanism, either control rod movement or turbine load 
change. These perturbations occur at fixed rates which limit the amount 
of reactivity that can be inserted during the simulation. 

Originally it was envisioned that an overall program would feature 
not only interactive input capability but also produce real time graphi- 
cal displays of the transient for a user at a time sharing computer 
terminal which had graphics capability. However, the computer language 
(IBM CSMP-III) required to solve the model's algorithms was not avail- 
able on the time sharing system, and furthermore, the programs long 
execution time precluded any real time response. 

In order to still provide as much user facility as possible with 
these restrictions, the interactive feature was partially retained for 
data input. Using a computer language available on the time sharing 


system (FORTRAN-IV), a program which interactively prompts the user to 
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enter specified controlling parameters was developed. This program in- 
corporates logic routines which permit inputs only within requested 
ranges. This feature allows for data reentry if a user error is made, 
and ensures that only inputs compatible with the model development's 
assumptions are used in the simulation. 

After the input is completed, the initial transient conditions are 
displayed at the terminal. A separate internal control program then 
transfers the user supplied inputs into the simulation program. The 
control program then transfers the simulation program to the batch 
processing system for execution. The interactive program notifies the 
user that this has been done. Hard copy plots of the transient are 
Subsequently produced. The post-shutdown Xenon behavior module is 
located in the interactive program and produces real time graphical 
displays at the terminal. 

Thus there are three distinct programs: 


1. An interactive program to prompt and receive user input, 
and also simulate post-shutdown Xenon behavior. 


2. A batch processed program which simulates the transient 
response to the user inputs... 


3. An internal control program which interfaces these two 
programs. 
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III. MODEL DEVELOPMENT 


A. POINT REACTOR KINETICS EQUATIONS 


A lengthly and formal derivation of the point kinetics equations is 
found in Reference 4 and will not be repeated here. More generally these 
equations are obtained from one group neutron diffusion theory with the 
assumption that the neutron flux is deparable in time and space, and 


with the inclusion of delayed neutrons [5]. They are listed below. 





d@(t) _ At)-A A 
at) ee) 2. 0,C,(t) + Q(t) (1) 
me. t ) fa 
j _ Sea - —_ 
at aK P(t) A.C. (t) Hinrse eae ote se (2) 
where @(t) = Flux amplitude function 
A(t) = Reactivity 
C, t) = Effective concentration of the ith delayed 
neutron precursor group 
Q(t) = Extraneous delayed neutron source strength 
ji = Effective delayed neutron fraction from 
the ith group 
A; = Decay constant of the ith group 
and > - 2.7; = Total effective delayed neutron 
7 fraction 


These equations are referred to as the point reactor kinetics equa- 
tions not because the reactor is considered as a single point in their 
application, but since spatial variations are neglected. f,, A, and A; 
are assumed to be constant. 

If the neutron flux spatial variation is assumed to be time invari- 


ant, Q(t) can be considered to represent the number of neutrons in the 
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core. Equation (1) may then be seen as a neutron rate equation where 
the three terms on the right hand side represent the rate of production 
of prompt, delayed, and source neutrons respectively. 

Equation (2) is a rate equation for the ith delayed neutron precursor 
group. The two terms on the right hand side represent production and 
decay rates respectively, The precursors are comprised of approximately 
thirty isotopes which historically have been divided into six groups 
with decay constants ranging from 0.0124 to 3.0 seconds a oe UN in 


the model, the precursors were considered to be represented by one 


effective group characterized by an averaged decay constant [2]. 


@ = 

(3: 

| 4 
[4 = = | 


=i 


and an effective delayed neutron fraction 


Ms 


f= ZA, 


u 


In equations (1) and (2), reactor power may be substituted for neutron 


flux if the precursor concentration is modified by C = E- 2 fold 


Since 
DeaEr eee 
where P- = Reactor power 
Ee = Energy released per fission 


Macroscopic fission cross section 


3M 
—h 
iN iN 


Integrated one group neutron flux 


For a reactor operating at power, the source term contribution to 


the overall] neutron population is negligible, i.e., Q=Q. 
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Incorporating these assumptions, equations (1) and (2) reduce to: 








pit) 2 Pave P(t) + Ac(t) (3) 
dc(t) _ 4 
fee =e P(t) - C(t) (4) 


1. Zero Power Reactor Transfer Function 
Consider an initially critical steady state reactor at some power 


level e at t $ 0. Equations (3) and (4) each yield: 


{3 
oo Oe (5) 
Now Tet P(t) = Br + o'P(t) 
C(t) =C, + dC(t) (6) 


A(t) =A, + Gt) 


where the zero subscript denotes the initial steady state value 
and the delta prefix a small perturbation imposed at t = 0 about 
this value. 

Noting that A, = 0 for a critical reactor and substituting equa- 


tions (5) and (6) into equations (3) and (4) yields for t 20 


d'P(t) = =. orn ee £. P(t) + AdgC(t) (7) 
& Act) = S(t) - act(t) (8) 


where the 2 term in equation (7) has been neglected. 


Taking the Laplace transforms of equations (7) and (8): 


s dP(s) = 2. ysis) = ZS PIS) + Ad'C(s) (9) 
/ 3 
C(s) = a FPS) sa Ad cts) (10) 


le 





Solving equation (10) for C(s), substituting this into equation 


(9), and rearranging, yields: 


d’P(s} zs R (11) 
a (S) S(A+ = 


This is the zero power reactor transfer function, so called because 
the reactor is assumed to be operating at a sufficiently low enough 
power that no feedback effects are realized. These effects are examined 
in the following sections. 

Equation (11) may be represented in block diagram form as shown in 


Figure 2. 


Jpis) SP is) 


Figure 2. Block Diagram of the Zero Power Reactor 
Transfer Function 


1 
fa 
(a+) 
(A + StxX 
2. Limitations of the Point Reactor Kinetics Equations 


The point reactor kinetics equations derivation is based on 


where Z(s) = 


the assumption that the spatial dependence of the neutron flux is 

negligible. This assumption limits the validity of these equations to 
transients where this remains a reasonable approximation such as those 
which result in only small changes in reactivity. The small perturba- 


tion technique used in the development of the zero power reactor 
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transfer function also requires only small changes in reactivity if the 
first order approximation is to hold. Specifically, must be less 
than 0.5/4. Physically, when yf is greater than 3 EMemnedceorr 115 
critical on prompt neutons alone. The simulated transients imposed on 
the model were limited to ensure that excessive amounts of reactivity 


were not introduced, 


B. REACTIVITY FEEDBACK MECHANISMS 


In the zero power reactor point reactor, the power level is assumed 
to be so low that it does not affect the reactivity, thus there are no 
feedback effects. However, for a reactor operating at a useful power 
level, feedback effects do exist and the reactivity becomes an implicit 
function of the reactor power level (or neutron flux). This dependence 
arises since reactivity depends on macroscopic cross sections which in- 


volve the atomic number densities of material in the core: 


> = NOs 


where 2s = Macroscopic cross section (cm~? ) 
N = Atomic number density (atoms/cm? ) 
GO = Microscopic cross section (cm ® ) 


The atomic number density can depend upon the reactor power level 
Since the concentrations of certain nuclei are constantly changing due 
to neutron interactions. Material densities also depend upon tempera- 


ture which is a function of reactor power level and hence the flux.’ 


iinderstadt. J. J. and Hamilton, L. J., Nuclear Reactor Analysis, 
Ist ed., Wiley, 1976. 


Ie 





| 


The two feedback mechanisms considered in the model were Xenon-let and 
moderator temperature. 

Reactivity can also be changed directly by an external source such 
as control rods containing a neutron absorbing material. Thus the over- 


all reactivity can be written as: 


= A564 Ax + Sr (12) 


where Ox, = Bee 
Se = FP 
Jd, = Overall core reactivity 
ide = Externally added reactivity 
fx = Xenon-135 feedback reactivity 
7dr = Moderator temperature feedback reactivity 
Fx = Xeonon transfer function 


Fy = Moderator temperature transfer function 
Equation (12) with the previously developed zero power reactor 
transfer function, given by equation (11), are incorporated in the 


block diagram below: 





Figure 3. Block Diagram of the Reactor Transfer Function 
with Xenon-135 and Moderator Temperature 
Feedback Loops, and External Reactivity 
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The next step in the model development was to derive the Xenon-135 
and moderator temperature transfer functions. 
1. Xenon-135 Transfer Function 
Xenon-135 is the most significant fission product poidon because 
of its enormous thermal neutron absorption cross section and relatively 
large fission yield. This isotope is not only produced directly from 
fission but also from the decay of other fission products as shown in 


Figure 4, 





Figure 4, Xe-135 Decay Scheme 


(1397) (135 


Since the / decay of Iodine-135 and Xenon-135 Xe), with 


the largest half-lives, are the controlling steps in this decay scheme, 


it was simplified by making the following assumptions: 


1) All ey is produced directly from fission (the production of 


Antimony-135 (13s pes 
instantaneous). 


Sb) and subsequent decay to I is considered 


2) The short lived metastable V35m 


3) The removal of i) by neutron absorption is negligible for 


Xe 1S ignored. 


the neuton flux levels used in the model (o!* neutrons/cm“sec). 
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iss 135 


feenl| the I decays to Xe. 


With these assumptions the effective decay scheme is shown in Figure 








oe 
Figure 5. Simplified Xe-135 Decay Scheme 
Using this decay scheme, the resulting rate equations for boo; and 
Pye are: 
Meee = yr Se Dee) ~ Az Tle) 
Teo = re Ot) —Axr (13) 
AX (é) 
Sa = Ya Lg Ble) + ATT le) Ax Xe) — H* Ble) X(t) (14) 
where lees 1997 number density (atoms/em>) 
X = 135y,, number density (atoms /cm*) 


ere. -1 
af = Macroscopic fission cross section (cm) 


Q = Integrated one group neutron flux (neutrons/cm“sec) 


ite Effective 135) fission yield 
yh Effective V39y, fission yield 
A, Sale| decay constant (sec!) 
melso -| 
Ay = Xe decay constant (sec) 
ics , 
a = Xe microscopic thermal neutron absorption cross 


section (cm2) 
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Again, using the first order perturbation technique, let 


ee 14 S(t) 
X(t) =X, + oAX(t) (15) 
Q(t) = 9, + Ap(t) 


where the zero subscript denotes an initial steady state value and 
the delta prefix denotes a small perturbation about this value. 


Upon substituting equations (15) into equations (13) and (14), and 


taking the Laplace transform, the relationship between the perturbation 


Mee and Pp is derived (See Apoendix A): 


SHS) _ Ba Bh-oa* Kos rz er BF +¥eEG-G" Xe) gig) (16) 
o's) S?7+ (os%Sotrx tAr)St Ar(GX Zot Ax) — 








The change in reactivity caused by a smal] perturbation in 139y, 
concentration is also derived in Appendix A: 
Phx (8) _ 4 
En oa are COX / (17) 
Xo + =U 
O3 
where Xo = Equilibrium USerye number density before the 
perturbation (atoms/cm) 
UY = Uranium-235 (2394) number density (atoms/cm°) 


oa Microscopic thermal neutron absorption cross 


section of oi) (en) 


A perturbation in neuton flux ( of’ )) is directly proportional to a 


perturbation in reactor power ( o’P) as shown by: 
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where Ee = Energy released per fission 








Za = Macroscopic fission cross section of 2391) (om~!) 
or 
SOs) _ 4 
= = k (18) 
PP(s) Eg dic 


The product of equations (17), (16), and (18) yields: 


ox x(S) of X(S) SG (Ss) Ax. (5) : 
dX (S) ASG (Ss) oP (Ss) oS P(s) x &x (S) K x tS ( ) 


Equation (19) is the transfer function relating reactor power and 
Fy. feedback reactivity. This equation is shown in block diagram 
form in Figure 6 where it has been incorporated with the previously de- 


rived zero power transfer function given by equation (11). 


S45) 


ape a pe 
+ 


Figure 6. Block Diagram of the Reactor Transfer 
Function with Xenon-135 Feedback and 
External Reactivity 


2. Moderator Temperature Feedback 


Reactivity feedback from changes in moderator temperature occurs 


as a result of changes with moderator density. The density is also a 





function of pressure, however, the pressure coefficient of reactivity is 
typically two orders of magnitude smaller than the temperature coeffi- 
cient, and was therefore not considered in the model development. The 
moderator density affects the moderator number density and hence the 
macroscopic scattering cross section of the moderator. 

The primary mechanism for the thermalization of the prompt and de- 
layed neutrons is by elastic scattering interactions with the moderator 
nuclei. Hence, variations in the macroscopic scattering cross section 
will affect the rate at which neutrons become thermalized. Changes in 
the thermalization rate affect the fission rate, or reactor power level. 
Power level changes affect the moderator temperature, thus a feedback 
loop is created. 

Because the moderator is also the coolant in the model, its temper- 
ature is not only a function of reactor power but also a function of the 
heat transfer process occurring in the heat exchanger, thus complicating 
the feedback loop. 

Therefore, in order to develop this feedback mechanism analytically, 
it is first necessary to model the plant's heat transfer processes. 

This thermal analysis is done in the following section. 

Another temperature feedback mechanism is the broadening of the 
Uranium-238 resonance absorption cross section for neutrons with increas- 
ing temperature. Because a highly enriched core was assumed in the 


model, with little Uranium-238, this effect was not considered. 
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C. THERMAL ANALYSIS 


1. Reactor Heat Transfer Function 
A lumped parameter model was assumed. This simplification pro- 

vided a set of ordinary differential equation which were sufficiently 
accurate for the simulated normal operating transients. In the lumped 
parameter model, heat transfer in the reactor was assumed to occur at a 
Single point. Thus, spatial variations were neglected. The core was 
considered to be a homogenized mixture of the uranium alloy fuel, the 
fuel cladding, and other structural materials, with a constant thermal 
capacity. 

The equation for the heat flow from the core was obtained from a 
basic heat balance. The heat generated from fission equals the heat 
required to change the temperature of the core materials plus the heat 


transferred to the coolant. On a per unit time basis, the heat balance 


1S 
dT, (t) 
ey = Ce may + hey | Te(t)- Tingy lB) | (20) 
where P(t) = Total Power generated in the core (Btu/sec) 
TA(t)= Average temperature of the core materials (°F) 


Tay(t) = Average reactor coolant temperature (°F) 


Oo 
tk 


F Total thermal capacity of the core materials (Btu/°F) 


h = Total heat transfer coefficient (Btu/° F sec) 


A similar heat transfer balance must also hold for the heat being 
transferred to the coolant and transported out of the core. This heat 


is the last term of equation (20) and is transferred to the coolant as: 
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hen | Tele) - Tavit) | = Cy Taye) + TimC [ Ta, (e) ~Te, te» | (21) 


where Cy = Total thermal capacity of coolant in the core 
(Btu/°F) 
ms = Coolant mass flow rate (1bm/sec) 
T, = Average reactor coolant outlet temperature cu 
0 
ie = Average reactor coolant inlet temperature (°F) 
C = Specific heat of reactor coolant (Btu/1bm°F) 


For simplification, the average reactor coolant temperature is 


assumed to be given by: 


Tav (t) = [ Tee) + Tite (&) | /2 B22 ) 


Rearranging equation (20) yields: 


oe d Trt) - T 74 ( 
Te (te) ay 14 aie = AV (x) +S P(e) (2s) 
~ Ce 
{ Pama ad 
where 1 in 


Substituting equation (22) into equation (21) to eliminate Th VeidGe 
8 





os — 
Iz a dlav(t) _ az 
TAY (t) [ire is | tle = regdia Tr(e} + 2 Te Ie, (t) (24) 
where To = ous ) and 
Myl 
Te= Su 
hem 


Zh 





Now let 


Te(t)= 7 ‘ + AT.(t) 
ct) = Tav.* Ste (25) 
Ty (ts Tho, oie 
Tp, (ts Tes, ae 


where the zero subscript denotes a steady state value and the 
prefix a small perturbation about this value. 

Substituting equations (25) into equations (22), (23), and (24), 
eliminating the average temperatures TA(t) and Tay(t) and taking Laplace 
Means forns, the following relationship between reactor coolant inlet and 


outlet temperature perturbations is derived (see Appendix B): 





72 7i 
-{ ef B(Z4)-7 1} s 1} Mey + Io P(s) 
STuo(s) = (26) 
7; Te fo (Zi 
{WB o5[ (Ze )+%]s+2] 
where 


, 








Cri2 MC 


This equation is the reactor heat transfer function for the reactor 
coolant outlet temperature as a function of the reactor coolant inlet 
temperature and reactor power. 

A block diagram representation of this transfer function is shown 


in Figure 7. 
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Figure 7. Block Diagram of the Reactor Heat 
Transfer Function 


In Figure 7 
_ (4%, [2 (Z,,)_7 
G.(S) 7 ( 2 S° + |Z (2 +1) -3]5+1] 
1 
G,(S) a 


2. Heat Exchanger Transfer Function 


As in the derivation of the reactor heat transfer function, a 


lumped parameter model was assumed for the heat exchanger. Two points 


of energy storage were assumed, the primary coolant water and the water 


On the secondary side of the heat exchanger. The thermal capacity of 


the heat exchanger metal was included with that for the secondary water, 


Since the thermal resistance on the primary side is predominant. A 


further assumption was made that the time spent by the primary coolant 
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while it passed through the heat exchanger was negligible in comparison 
with the time spent in the primary piping. 

The state of the steam produced on the secondary side of the heat 
exchanger was assumed to always be dry and saturated and that the 
secondary water is always at the saturation temperature for the existing 
pressure. These assumptions were justified because moisture separators 
and recirculation can provide high quality steam and preheating of the 
feedwater. 

With these assumptions, the following equations were obtained from 


a heat balance per unit time: 


; TV 

Mm | Zi, (t) - Te, (a) | = Cm vets) + him | Tavg(t) - Ts (¢) | 6273) 
hm | Tavg(e -Tse) ] = Cs SBE + RG) (28) 
where T(t) = Heat exchanger coolant inlet temperature (°F) 


Teo(t) = Heat exchanger coolant outlet temperature (°F) 


Tay (t) = Heat exchanger average coolant temperature (°F) 
B 
T.(t) = Saturated steam temperature (°F) 
Ney = Total heat transfer coefficient (BTU/°F.sec) 
Cu = Total thermal capacity of coolant in the heat 
exchanger (BTU/°F) 
C. = Total thermal capacity of heat exchanger metal 


and secondary water and steam (Btu/°F) 
P(t) = Power delivered by the heat exchanger (BTU/sec) 


m = Coolant mass flow rate (1bm/sec) 
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The total heat transfer coefficient, a function of the primary coolant 
flow rate (a constant in the model) and the heat transfer characteristics 
of the heat exchanger, was assumed to be constant. Also, the time delay 
in transferring heat across the heat exchanger tubes was neglected. 

This changed the shape of the initial thermal transient but had little 
effect on the basic dynamics of the secondary loop. 

Again, as a simplification, the average temperature of the coolant 


in the heat exchanger was assumed to be given by 


— 
pig0c a2 


(29) 
The power delivered by the steam generator is proportional to the 


product of the steam flow rate and the difference in enthalpy between 


the steam and feedwater 





hm [Tava (&) - Ts)| = mg (Hs-Hew) + Cs ott 


where mM. = Steam flow rate (1bm/sec) 
ie = Saturated steam enthalpy (Btu/1bm) 
He = Feedwater enthalpy (Btu/1bm) 


This equation assumes the steam and feedwater flow rates are always 
equal and thus neglects any instabilities in the secondary steam. 

As discussed in Reference 4, the enthalpy of saturated steam is 
nearly constant over a wide range of pressure, varying from 1198 Btu/1bm 
to 1204 Btu/lbm over a pressure range from 200 to 800 psig. The feed- 
water enthalpy,which depends on condenser pressure, is usually between 


50 and 100 Btu/lbm. Thus the enthalpy difference may be regarded as a 
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constant, and the power delivered by the heat exchanger is directly 
proportional to the steam flow rate. 

The impedance to steam flow caused by the turbine is nearly indepen- 
dent of turbine speed. If constant backpressure is assumed, the steam 
flow rate is directly proportional to the throttle opening at a given 


pressure. Thus, 


mM. = gah 
where TF . = Saturated steam pressure (1bf/in®) 
A = Proportionality factor which is_a function 


of the throttle setting Ibm ine 
bt sec 


The numerous assumptions made in the heat exchanger model development 
limit the accuracy of the resulting equations. However, the errors in- 
volved in these assumptions are usually less than the amount of uncertainty 
in the engineering valves of the coefficients used in the equations. 

Recalling the previous assumption that 

dT ay (ice 
eee = -() and 
dt 


Substituting equation (29) into equation (27), yields: 


mam [Thy (t) — Teo(e) | = him { | Tay (e + Teo) | - Tste) | (30) 


Solving equation (30) for Te , gives: 
0 


recent. (1k 
Om oe ee Wy (2-Ka ) (31) 
(1 +k) 
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2mmC 


where K, = 
, 1 hm 


Substituting equation (29) into equation (28) gives: 


~— hm F ocr Tete) | 500 | = Cs 2s + P(t) 


Rearranging this expression, 


d Ts (t) 
- — + Tse) = 4 | Ta, ce - Teo e| + = Fi(e) (32) 
Let amet + oT, (t) 
10 1 
ite Ue) Sune Aue, (33) 
0 05 
Te(t) = Ty + JT.(t) 
ane ae eee) 


where the lowest zero subscript denotes a steady state value and 
the delta prefix a small perturbation about this value. 

By substituting equations (33) into equations (31) and (32), and 
taking the Laplace transform, the following set of equations is derived 


(see Appendix C): 


> | Tags) +o'Tus (3) | at AP, (S) 





Te = 34 
d'Ts (S) ar (34) 
yen 20°Ts (S) — o Tus (5) ce Ka | (35) 


1+ Ki 


acs 
where ie = cara 





These equations are represented in a block diagram in Figure 8. 


dT (s) 





ro Ps) 


Figure 8. Block Diagram of the Heat Exchanger 
Transfer Function 
With oT, and oP, as inputs, equations (34) and (35) form an 
1 
algebraic loop in oar and 71: These two equations were solved 
0 
with the use of an inherent functional routine available in the CSMP-III 


language which was used to formulate the model. 
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3. Primary Piping Transfer Functions 


While circulating through the primary loop, the coolant under- 
goes mixing and transport delay effects. For example, a transient in 
the coolant temperature at the heat exchanger outlet does not appear at 
the reactor inlet until sometime later. Where there are volume or flow 
direction changes, as in the reactor coolant inlet plenum, mixing occurs 
causing a smearing of a temperature transient. Both of these effects 
were approximated in the model by combinations of two types of time 
delays: a pure transport delay and a simple mixing delay. 

Assuming no mixing in the primary piping and no heat loss with 
perfectly insulated pipes, pure transport delays are encountered in the 
Diping runs between the reactor and heat exchanger. With these assump- 
tions, the temperatures involved in the transfer of heat between the re- 


actor and heat exchanger can be given as 


(t a T,) 


where the inlet coolant temperature to the heat exchanger ey has the 
1 


Same form as the outlet temperature of the reactor Ty after a fixed 
0 
transport delay T3. This transport delay can be approximated by the 


following differential equation which is derived in Appendix D. 


Ine th Tes i = T 


4 (36 ) 


j dt No 


Similarly, the transport delay from the heat exchanger outlet to 


the reactor inlet plenum is: 


(t) of alm qT (t) (37) 


T 2 
25 ip = Ty (t) 


C 





where Te. (t) = Reactor inlet plenum coolant temperature (°F) 


4 
© 
~-——s 
ct 
ee 
it 


Heat exchanger outlet coolant temperature (°F) 


This can be approximated by 


pe dTe., (t) 


eo + |4 pane = Jog(t) (37) 


As in the previous development, let 


H 
09 0 
Wy Ve etre ae ale ee (38) 
1p 1po 1p 
fet) = it. + ot. (t) 
Se Oh oe 
eC lipueat | oT, (t) 
1p 1po Ip 


where again the lowest zero subscript denotes a steady state value 
and the delta prefix a small perturbation about this value. 

After substituting equations (38) into equations (36) and (37), 
taking the Laplace transform, the following transfer functions are 


derived (see in Appendix D). 


d'Tiisp (S) sl 
—<—<—— (39) 
d’Tu, (s) 1+T35 


d'Teg(S) 1+ 745 
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Combined mixing and transport effects are encountered at the reactor 
and heat exchanger inlet and outlet coolant plenums. For simplification, 
only the mixing effects at the inlet plenums were considered in the 
model. Assuming perfect mixing and after performing a heat balance on 
the plenum concerned, the combined mixing and transport effects are ex- 


pressed by differential equations developed in Appendix D of the form 


7 HME) 4 T(t) = T(t) 


0 
where 
T.(t) = Plenum inlet temperature (°F) 
TA ft) = Plenum outlet temperature (°F) 
7 = Mixing delay 


As before, by using a small perturbation technique, and taking the 
Laplace transform, the following transfer function is derived (see 


Appendix D): 


d Tio (S) Sey a (41) 
PTH: (Ss) 1+7S 


Block diagram representations of the coolant piping transport and 


mixing transfer functions is shown in Figure 9. 


D. OVERALL PLANT BLOCK DIAGRAM 


The transfer functions developed previously for the zero power point 
meactor, re feedback, reactor and heat exchanger heat transfer, and 
the primary piping were interconnected resulting in the overall model 


block diagram shown in Figure 10. The moderator temperature feedback 


ay 





Jd P(s) 





IT ..(5) 
REACTOR 






HEAT 


STAs) | EXCHANGER 


PAs) 


Figure 9. Block Diagram of the Transport 
and Mixing Delay Functions 


loop is seen to consist of the heat transfer process and primary piping 


transfer functions which yield Za and 7, . hese temperatures 
i 0 


are summed, then halved, yielding Taye which is then multiplied by 


the negative temperature coefficient ae generating the moderator 


temperature feedback reactivity jo 


MODEL CONSTANTS 


The following plant parameters from the Shippingport Atomic Power 


Station were obtained from Reference 6 and used in the model development. 
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Overall Plant Block Diagram 
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General Parameters 


Reactor thermal power 231 MW 
Reactor coolant system oressure 2000 psia 
Reactor coolant average temperature Sao 
Steam pressure at full load 600 psia 


Reactor Coolant System 


Reactor coolant flow rate 6280 1bm/sec 
Reactor coolant outlet temperature Geo er 
Reactor coolant inlet temperature B0o) 
Coolant volume in core 103 ft? 


Reactor Core 


Configuration Right cylinder 
Size Crom pee Gide sO flees ming 
Fuel load uc? (seed) oeikG 


Composition (seed) 


1) Water 32500 
2) Fuel alloy 30.0 v/o 
Sve iincal toy, Cae 


Control Rods 
1) Total rod worth 0.256 k 


2) Scram time 0.35 sec time delay 
1.0 sec rod drop 


Nuclear Data 


dad. 


D’: 


Thermal neutron flux Zak 10/4 nfem® SeC 


, =5 
Prompt neutron lifetime 5.6 X 10 ~ sec 
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c. Effective delayed neutron fraction 0.0077 


d. Temperature coefficient of reactivity Olek 1077 SP k/°oF 


5. Reactor Protection Setpoints 


a. Scram 
1) High reactor power 138% 
2) High reactor coolant outlet temperature 550°F 
b. Cutback 
1) High reactor power 114% 
2) High Startup rate 1.74 Decades per min. 
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A. INTERACTIVE PROGRAM 


Figures 10-16 are examples of the interactive program output. As 
Shown in Figure 10, the user is first given a description of the purpose 
of the overall program and general instructions for input entry. As 
Shown in Figure 11, the user is then prompted to enter an initial steady 
State power level, after which the corresponding major plant parameters 
are displayed. The user is then prompted to choose the type of simula- 
tion, either plant transient at power, or post-shutdown Xenon-135 
behavior. The user must be at the Tektronix 4012 console in the computer 
center if the latter is chosen. 

If the post-shutdown Xenon-135 behavior is chosen to be examined, 
graphs such as those shown in Figure 13 are generated and displayed on 
the Tektronix's screen. In addition to showing the time response of the 
Xenon, the time of its peak and associated maximum reactivity are dis- 
played. 

If the plant transient simulation is chosen, the user is then prompted 
to choose the mechanism for initiating the transient, either control rod 
Movement or turbine load change. In Figure 14, the user has chosen the 
former. The program then prompts the user to enter the direction and 
time length of the control rod movement. In Figure 15, the turbine load 
change was chosen and the program request the final turbine load. In 
either case, a summary of the transient inputs is displayed with a 


notice that the interactive portion is complete. 
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Figure 13. Post-shutdown Xenon-135 Behavior 
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Figure 16 is an example of the interactive program's logic routines 
which will only permit inputs within the requested ranges to be accepted 


for the simulation. 
B. SIMULATION PROGRAM 


The program generates Versatec plots showing the plant's power, re- 
activity, and temperature transient responses for the user's inputs. 
Figures 17-20 are composites of these plots. 

In Figure 17, the transient was initiated by a simulated 20 to 40 
percent ramp change in turbine load at 1/2 percent per second. In 
Figure 18, the transient was initiated by a simulated 60 to 40 percent 
Change in turbine load in the same manner. Both of these figures show 
the characteristic power demand following and inherent stability response 
typical of a PWR plant with a constant average temperature program and 
negative temperature coefficient. 

The inherent stability feature is further displayed in Figures 19 
and 20. In Figure 19, the transient is initiated by a simulated 10 second 
inward movement of the control rods at a reactivity insertion rate of 
mic X ery Sk per second from an initial 50 percent power level. In 
Figure 20, the transient response to an outward control rod movement 
with the same parameters is shown. As seen in both cases, the reactor 
Dower stabilizes about the initial turbine load and the average reactor 
coolant temperature stabilizes at a new level to compensate, via the 
negative temperature coefficient, for the reactivity inserted by the 


eontro! rods. 
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Figure 17. Power, Reactivity, and Temperature Response 
to a 20 to 40% Turbine Load Change 
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Figure 18. Power, Reactivity, and Temperature Response 
to a 60 to 40% Turbine Load Change 
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Figure 19. Power, Reactivity, and Temperature Response 
to a 10 second Inward Control Rod Movement 
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Figure 20. Power, Reactivity, and Temperature Response 
to a 10 second Outward Control Rod Movement 
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V. CONCLUSIONS AND RECOMMENDATIONS 


The model responds to normal operating transients with the charac- 
teristics of a PWR plant. The interactive program provides the user with 
the facility to initiate the simulation or examine a real time display 
of post-shutdown Xenon-135 buildup and decay. 

As previously discussed, the purpose of this work was to develop a 
learning device for nuclear engineering students at the Naval Postgrad- 
uate School. The result was a relatively simple model of a PWR power 
plant and a simulation program that suffers from a long execution time. 

The model's simplicity resulted mainly from the consideration of 
an effective single group of delayed neutron precursors in the reactor 
kinetic equations and the use of lumped parameter analysis in the plant's 
heat transfer processes. The model's sophistication could be increased 
by: 


1. Expanding the reactor kinetic equations to consider six 
delayed neutron precursor groups. 


2. Developing a multi-section heat transfer model for the 
reactor and the heat exchanger. 


3. Incorporating variable reactor coolant mass flow rate and 
heat transfer coefficients. 


These additional features will also increase the complexity of the 
simulation program, thus aggravating the long run time problem. This 
concern might be eased by taking a different approach to the model's 
formulation than the transfer function method. While this method was 
readily coded using the CSMP-III language, physical insight to the plant 
dynamic processes was lost when the Laplace transform and subsequent 


grouping of constant terms was performed. The use of state variable 


54 





theory might not only allow this physical appreciation to be retained 
but also result in shorter program run times. 

Regardless of the formulation method, a real time response to the 
Simulation can probably only be achieved by using an analog computer. 
The wide range of time constants associated with the equations (prompt 
neutron life-times on the order of 107° seconds to Xenon decay half 
lives on the order of hours) require a small numerical integration time 
interval over a large time period. The result is a prohibitively long 
run time on a digital computer for an interactive program. 

It is hoped that refinement of the model will be continued. While 
the existing model does reflect the general characteristics of a PWR 
plant to normal operating transients, it does exhibit relatively large 
power over/under shoots in response to turbine load changes as seen in 


Figures 1/7 and 18. 
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APPENDIX A 


DEVELOPMENT OF THE XENON FEEDBACK TRANSFER FUNCTION COMPONENTS 


1. Derivation of G(S) 


The previously developed rate equations for 285 and eae from the 
effective Xe decay scheme shown in Fiqure 5 are: 

sf = Sr UeSlt) — Arle) (A 1) 
ak = 3x 2.£DB (te) + ArT (te) —AxX(e) — 09% Dee) X(t) (A 2) 
dt 
where I(t) = Vo97 number density (atoms/cm?) 

R(t) = eee number density (atoms/cm?) 

2.f = Macroscopic fission cross section of 7 Com”! ) 

Q = Average integrated one group flux (neutrons /cm“sec) 

3; = Effective beey fission yield 

Sy = Effective! Xe fission yield 

AT = eT decay constant (sec™!) 

Ay = '3%y¢ decay constant (sec!) 

om = oe thermal neutron absorption cross section (cm?) 
fee I(t) = 1, + oi(t) 

MEEK, + ox) (A 3) 


Q(t) = 9 + AD(t) 


where the zero subscript denotes an initial steady state value and 


the delta prefix denotes a small perturbation about this value. 
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an - i 


7 i 


With the reactor at an initial steady state with a flux level on the 


is SS 


equilibrium values of I and Xe are found by setting their time 


dependence in equations (A 1) and (A 2) equal to zero. Thus, 


yi 
= Jt zh®o (A 4) 
AL 
Xo = ¥«x ZF Bo + Arlo (A 5) 
AMX + SeaGo 
Substituting equation (A 4) into equation (A 5) 
x, = (x4 32) BE Bo (A 6) 


Ax + 3 Bo 


Substituting equations (A 3), (A 4), and (A 6) into equation (A 1) 


and (A 2), neglecting the dg dX term (first order approximation) 


il 


Sot Yr Be /D-AroL (A 7) 


SX = AgoT + (3x Ze - Xo) SD - (Ay +O Do) AX (A) 


a 
St 


Taking the Laplace transform of equations (A 7) and (A 8) and solving 
for oI(S) and dXx(S) 


¥x Ze /SDS is) 
fI(s) = a (A 9) 
= Ga.Xa 
iii. Ard'I(s) + (¥x ZP-o ) SSos) (A10) 


CStAx + O5'Zo) 


Si, 





Substituting equation (A 9) into equation (A110) 
At Ix ZeSSHGS) + (star) (¥x Le - 53% Xo) SDs) 





SX(s) = 
CS5#Ar) (StAx+ 05*%Gq) 
dx(s) 
Expanding, collecting terms, and solving for Gigs) 
SKS) = (3x EE ~O3 Xo) St Ar (¥r Se + Ix LE -g* Xo) ene 


This is equation (16) on page 23. 


2. Derivation of x 
The effective core neutron multiplication factor k is given by 


the familiar “six factor formula" found in the literature 


k “Rept PeNL PINE eh) 
where 4 = Thermal fission factor of the fule 
€ = Fast fission factor of the fuel 
p = Resonance escape probability 
f = Thermal utilization factor in the core 
PENL = Fast neutron non-leakage probability 
Pov = Thermal neutron non-leakage probability 


The core reactivity is defined as 
_ kel 
Go (A 12) 
Consider a reactor at an initial critical steady state condition 


ss 


with an equilibrium Xe concentration Xo. By the definition of 


eeiticality 





ko 
Po 


Heol 
aS 


(13) 
In this condition 


f,= Ze. __%3 


a a 


core SF, Xo (A 14) 
Bere Pee = Fuel macroscopic thermal neutron absorption cross 
section 
pa = Initial core macroscopic thermal neutron absorption 
do cross section 
Xo _ _ 35 : 
ee = Initial Xe macroscopic thermal neutron absorption 


cross section 


and the macroscopic thermal neutron absorption cross sections of the 
moderator and other core materials has been neglected, a reasonable ap- 


proximation in a highly enriched core. 


eS 


Now impose a small perturbation in the Xe concentration on the 


core. Neglecting the effect of the Ise 


a reasonable assumption for a large reactor, 


X=Xo + dX 
k=ko+ ok 
f = fo + df 


P=fPo+ IP = AO since = 0 
The remainder of the terms in equation (A 11) are unchanged, thus 


dk = of 


and 


ve k since ko = 1 
fo 


ALA 
oO 


99 


Xe perturbation on neutron leakage, 





where 


F F 
se fs ee pen) 
5 2g, + Oa 0X 


from equation (A 12) and (A 13) 
di = Sy = e-1 | 4 mde (A 16) 


substituting equations (A 14) and (A 15) into equation (A 16) 





Sox = mS poso'x Becerra 
Dae er Cae pe 
or 
JPx = = 4 = X% 
ox Xot 23. Su 


This is equation (17) on page 23. Since fuel depletion effects are 
not considered during the short simulated transients, the fuel number 
density U is assumed constant and Sy is invariant in the time and the 


Laplace domains. 
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APPENDIX B 


DEVELOPMENT OF THE REACTOR HEAT TRANSFER FUNCTIONS 


The previously developed equations describing heat transfer from the 


reactor to the coolant are 


+ = 
Tecey + Ta SEE = Taye) + SE PEE) (B 1) 
if (e) rf 
Tav (ey) 2 + 2 = +t T2 cnt = Tr(t)+2 es (t) (B 2) 
Tav (4) = 3 (thet) + Tes () | (B 3) 
where Te(t) = Average core material temperature (°F) 


Average reactor coolant temperature (°F) 


“4 
> 
<= 

— 
ctr 
—S 
il 


Ty (t) = Reactor coolant outlet temperature (°F) 


T_.(t) = Reactor coolant inlet temperature (°F) 





Ci 
P(t) = Total power generated in the core (Btu/sec) 
Ce = Total thermal capacity of core materials (Btu/°F) 
Cu = Total thermal capacity of coolant in the core (Btu/°F) 
E = Specific heat of coolant (Btu/bim°F) 
m1 = Coolant mass flow rate (1bm/sec) 
New = Total heat transfer coefficient (Btu/sec°F) 
and 
Pe Cm 
lo = e 
mC 
hem 
Te = Sn. 
hem 
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Let P(t) =P. + oP(t) 


0 
a) meee eT (£) 
Tyy(f) = Tyy #0 AT ay (t) (B 4) 


10 
where the lowest Zero subscript denotes a steady state value and 
the delta prefix a small perturbation about this value. 
From substitution of equations (B 4) into equations (B 1), (B 2), 


and (B 3), with the reactor in a steady state at t $0, 
TF, = Tivo + 810 
-_ 
see j2 
Tavo (4 +2 fey = Tro + 2 4 Te;, (B 5) 
at 
Tavo = 4% (Tog + Tes.) 


and 
SP(a) =o 


J TE (o) = dD Tay (o) = 'Tu,(0) = ole; (oO) = 


Impose the perturbations at t=0. Substituting equations (B 4) and 


(B 6) 


(8 5) into equations (B 1), (B 2), and (B 3) with the initial conditions 


given in equations (8B 6), 


o'Te (te) + Tt ¥ PTe(e) = J’ Tay (€) eS ce HPC) (B 7) 


T: a T: 
d'Tav (t) [sez =| +72 $ | ota) | = Tet) +2 - oTex(t) — (B 8) 


oS’ Tav(t) > 5 | Tae (t) + OTe, ce) | (B 9) 
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Taking the Laplace transform of equations (B 7), (B 8), and (B 9) 


with the initial conditions given by equations (8B 6) 


Te (Ss) [ sT4+ 1 | = Tav(S) + ce PCs) (B 10) 
12 = Tz 

Tav(s) | ste+142 de | = Te(s) + 2 = Te (Ss) (B 11) 

Tavis) = 3 | Tis) + Te (3) (B 12) 


where the delta prefix and lowest subscript have been deleted for 
readability 


Substituting equation (B 12) into equations (B 10) and (B 11) 


Te (si4+1) = = (th + Te) Pas =P (Eas) 
~ 72 T: 
Or.) (site reef) » = Te +2 Te (B 14) 


where the s domain dependence notation has been deleted for 
readability. 


Solving equation (B 13) for Te and substituting into equation (B 14) 


—_ 
(Tut Te) (stot1+2 = Iz) ( Ta+1) => (Tu+Te) + - P+2 = (st4+ 1) Te 
After expanding the products and collecting terms 


Tr Tz, , it, Taiz T; i112 71 _ Tae 12], 7% 
Th =" 2 +( $+ 2 + Bie) 54% |= - Te [s* & peek. 2) s- Els Bp 





After multiplying through by and solving for Ti a (s)) 
0 


ee 

b. _ (BBs +[ E(ap+ 1)-Ta]s-a1 Pres) + ¥ PP (s) 
= 
{s+ [R(ge)ea]se| 
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a 
In lo 


—; and all notation has been restored. 
Sp 12 


where y= 


This is equation (26) on page 28. 
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APPENDIX C 
DEVELOPMENT OF THE HEAT EXCHANGER TRANSFER FUNCTION 


a 


The previously derived heat exchanger heat transfer equations are 


o Tete) - THs (+) [4- Ka] 





: (C1) 
Teg (©) | 4 + K+ | 

1 
Ss diste) tTs(t) =F | Thay cer + Teo tt) | + Ary Putt) G2) 
hm™ dt 
where Ty, (t) = Heat exchanger coolant inlet (°F) 


Te, (t) = Heat exchanger coolant outlet (°F) 





T.(t) = Saturated steam pressure (PSIQ) 
P, (t) = Power delivered by heat exchanger (Btu/sec) 
Co = Total thermal capacity of heat exchanger 
metal and secondary water and steam (Btu/°F) 
C = Specific heat of coolant (Btu/1bm°F) 
Nas = Total heat transfer coefficient (Btu/°F.sec) 
Mm = Coolant mass flow rate (1bm/sec) 
and 2myC 
he 
Let Twilt) = Ty Oeinen) 
Pome oT) (Gas) 
T(t) =T, + OT.(t) 
2 Co) Pi me 2) 


where the lowest zero subscript denotes a steady state value and the 


delta prefix a small perturbation about this value. 
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With the plant in an initial steady state at t <0, equations (C 1) 
and (C 2) are 


21s, m= THa, [4+ K1 \ 
ete One 


4c - | (C 4) 

Con [1+ K4 | 

iliceee = (Ths. + Teo.) a ~~ Hie (co) 
and 

d'F.4 (9) = © (C 6) 


J’ Tu: (0) = o’Te, (0) = oTs(o) = 
Impose the perturbations substituting equations (C 3), (C 4), and 


Memo) into equations (C 1) and (C 2) 


E. ad'Ts5 - dH; (1 -K12) 
i.  (Gsx:) (C 7) 
Cs d x 
a se 0 Ts (t) * oT Ce) = + | Tis (e) + JTeo(t) | _ — JR (e) Gre) 


Taking the Laplace transforms of equations (C 7) and (C 8) with the 
initial conditions given by equation (C 6) 
5 | PTeg(s) + ATH, (Ss) | - FE, SPLIS) 
2 eS a Hs (5 — W+m tS 


d Ts (Ss) = 
TIsSti 


where 75 = Csi 
hm 
7'Te,(8) = {20°TS(8)-o'TH, (s) [1-x1] } / [1+k4] 


These are equations (34) and (35) on page 33. 
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APPENDIX D 
DEVELOPMENT OF THE PRIMARY PIPING TRANSFER FUNCTIONS 


1. Transport Delay Transfer Functions 


The transport delay of the coolant between the outlet of the reactor 


and the heat exchanger inlet plenum can be expressed as 


TH) See emlia) (D1) 


= 
= 
m 
“ 
© 
of 
ate 
te 
——, 
ctr 
eee” 
i 


Heat exchanger inlet plenum coolant 
temperature (°F) 


Reactor coolant outlet temperature (°F) 


—i 
Xx 
Ww OO 
one, 
ctr 
eel 
i ib 


Transport time delay (sec) 


Rearranging terms and expanding equation (D 1) in a Taylor series 


- iSite) 
To ©) = Thy, (t) + 13 Sieh + — (eZ) 


For slow temperature changes, the second order and higher terms 


in equation (D 2) may be ignored. 
» dk 
Tay, ce) + 13 ae’ = Tho (4) 


Let 


wilt) = TH+ Pty, Ce) (D 4) 


ie dT, (t) 


eT) 


T 


Pja=< 


Ho! 


where the lowest zero subscript denotes a steady state value and 


the delta prefix a small perturbation about this value. 


6/7 





From substitution of equations (D 4) into equation (D 3) and for 


the temperatures in a steady state at tS0 


eT 
ites Ho, (Dies) 


and 


ATH) = STHo(t) = 0 us 


Impose the perturbations at t=0. Substituting equations (D 4) and 
(D 5) into equation (D 3) with the initial conditions given in equation 


(D 6) 


J Tus, (t) + 13 £- sure) = moma) (D 7) 


Taking the Laplace transform of equation (D 7) and solving for 


ee / JT... 





PTHyp(S) if 
DP TH9 (3) 1 + Tas 


This is equation (39) on page 36, 

By following a similar development, the transport delay between a 
perturbation in the heat exchanger coolant outlet temperature and the 
resulting perturbation in the reactor inlet plenum coolant temperature 


can be expressed by the transfer function 


dS Teip (s) ae 4 


J‘Tes (S) . 14+ 145 





where Plc; 5 = Reactor inlet plenum coolant temperature 
perturbation (°F) 
dS'Te g = Heat exchanger coolant outlet temperature 
perturbation (°F) 
a : 
1, = Transport time delay (sec) 
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This is equation (40) on page 


2. Mixing Delay Transfer Function 


Let the reactor and heat exchanger coolant inlet plenums be represent- 


ed by the control volume shown in Figure D-1. 


M = Mass of coolant in plenum (1bm/sec) 


m,To 

T = Average plenum coolant temperature 

ory) 

T. = Plenum inlet coolant temperature 
. ae, 

T. = Plenum outlet coolant temperature 
ee a(er) 

ro m = Coolant mass flow rate (1bm/sec) 


Figure D-1. Mixing Volume 


If perfect mixing is assumed, the coolant outlet temperature is equal 
to the average coolant temperature in the plenum. At heat balance on 
the volume requires that the net heat flowing into the plenum be equal 
to the increase in stored energy. Assuming no ambient heat losses, and 
assuming the temperature of the plenum structural material remains 


constant, the heat balance per unit time over a time interval At is 








d T(t) oJ To (t) 
mAtC [ Tce) - Tote) | = MC ~y-— St = mate ee 
where C = specific heat of coolant 
The resulting equation is 
mM ed To (t) » dIo(t) _T. 
ans = +Ta(t) = Tj tt) (D 8) 
a aS + To (te) J ae o ) 1 
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where T = — == mixing time delay. 


Let 


+ 
~~ 
ct 
— 
il 


aM) 
ice Som mca) 


where the lowest zero subscript denotes an initial steady state value 
and the delta prefix a small perturbation about this value. 


Substituting equations (D 8) into equation (D 9) 
7 ZS Toe) + Soe) = ST Ue) (D10) 
Taking the Laplace transform of equation (010) and solving for 
dT./ dT. 


dT, (s) 
fT. (s) 7 has 


T 


This is equation (41) on page 37. 
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APPENDIX E 
INTERACTIVE PROGRAM LISTING 
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